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Part I. Discussion of the Extended W.K.B. Approximation 

The W.K.B. approximation for the solution of the height-gain differential equation for 
a curved stratified troposphere is discussed in detail. The approximation depends mainly 
on a variable ui(r) which can be interpreted as the height dependent contribution of the 
phase for a field solution obtained by separation of variables. An expansion of in(r) with the 
aid of partial integrations leads to further approximations which facilitate the determination 
of the eigenvalues, and of the amplitudes of the modes connected with the propagation 
problem. The influence of the refractive-index profile, if assumed as smooth, then appears 
to be restricted to a dependence on the surface values of this index and of its gradient insofar 
as propagation over the ground is concerned. Further, all height effects of elevated antennas 
can be expressed in terms of the distance to the corresponding radio horizon. This results 
in simple relations between the fields connected with two different refractive-index profiles, 
provided both profiles coincide near the earth's surface. 

1. Introduction 

The propagation theory of a concentrically stratified atmosphere usually concerns a 
discussion of the corresponding height-gain differential equation. The W.K.B. approximation of 
the latter has amply been discussed, but little attention has been paid to its corrections. Some 
material on the form of these corrections has been presented by Pekeris [l] 2 while applying 
an earth-flattening approximation. This paper concerns (in part II) an expansion of the 
solution for a curved stratified atmosphere which starts with an "extended W.K.B. approxi- 
mation." By the latter we understand a well-known modification of the W.K.B. approxi- 
mation, in terms of Hankel functions of order 1/3; it is determined such as to remain finite at 
a turning point of the geometrical optics trajectory associated with the solution. The dis- 
cussion of this extended W.K.B. approximation (in part I) shows the dominating role of the 
refractive-index profile near the earth's surface, provided that this profile and all its derivatives 
are continuous functions of the height throughout the troposphere. 

2. Reduction to a Scalar Problem 

We start from Maxwell's equations for time-harmonic solutions (time factor e~ iuf ) for a me- 
dium with spherical symmetry, the refractive index n(r) of which only depends on the distance 
r to the center of symmetry. We further assume a permeability 1; the equations in question 
then read as follows in Gaussian cgs units: 

Curle-i-h=0, (1) 

c 

Cuvlh+i-n 2 (r)e=0. (2) 



i Contributions from Philips Research Laboratories, N. V. Philips' Gloeilampenfabrieken, Eindhoven, Netherlands. 
2 Figures in brackets indicate the literature references at the end of this paper. 
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The two standard solutions, known as the electric and the magnetic solution respectively, 
can be represented by: 

e e= ^r cirri curl {n(r)H e r}> 

w n\r) y (3) 

h e =curl {n(r)U e r}, 



"} 



and 

loo 



•curl {n 2 (r)U m r} 



^ cn\r) 

h m =curl l^^y-curl {^ 2 (r)n w r}J- 



(4) 



The symbol r represents the radial vector of length (x 2 +y 2 + z 2 ) 1/2 . 

The Maxwell eq (2) is verified at once for the electric solution (3); the same holds with 
respect to the other Maxwell eq (1) for the magnetic solution (4). The remaining Maxwell 
equations can be checked by deriving the following relations with the aid of a tedious analysis 
based on vector-field identities: 

Curl e e — - h e = curl r4 — AIT+(-^ — r? - )!! >> 

c o) [ ^ V dr z n/ J 



Curl h m +—n 2 e m 



- curl |r/ AII+(- 2 - n 2 U \~l 



Therefore, both Maxwell equations are satisfied, for the electric as well as the magnetic 
solution, if the scalar II does satisfy a wave equation of the form, 



An+ ? r&(r)II=0; (5) 



the effective refractive index then has to be defined as follows: 

c 2 d 2 1 
fteffM =n2 — 9 n--7-^ - for the electric solution. 
e K J co 2 dr 2 n ; 

( 6 ) 
nl s (r)=n 2 for the magnetic solution. 

3. Height-Gain Differential Equation and Its Extended W.K.B. Approximation 

We are particularly interested in the solutions corresponding to a vertical electric or magnetic 
dipole. In the system of spherical coordinates r, 6, p(x-{-iy=r sin 6e i<p , z=r cos 6) this dipole 
may be situated at r=b, 6=0. The field then becomes independent of <p, and particular solu- 
tions Hi of (5) are found by a separation of variables according to 

U l =f l (r)P l (cos d). (7) 

The corresponding height-gain differential equation for/ z (r) can be put in the form 

{^+^mf(r)\{r/ I (r)}=0 ) (8) 

if 

m?(r)=nyr)-&I), (9) 

k =u/c being the wave number in vacuum. 

The W.K.B. approximation of (8) consists of a linear combination of the functions 

p±ik C r mj(s) ds 
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n l herefore breaks down at a so-called turning point, that is a zero r—ri of mi(r). In the case 
of a single turning point a better approximation is obtained, as is well known, by approximat- 
ing m 2 ,{r) near r=ri by a linear profile; the equation can then be solved there rigorously with 
the aid of Hankel functions of order 1/3. This procedure leads to the following "extended 
W.K.B. approximation:" 



rf l (r)^A l 



\ I. ml(s)ds I f . 



it reduces to the conventional W.K.B. approximation if the Hankel function is replaced by 
the first term of its asymptotic expansion for large arguments. 

In (10) we assume Im ra z >0, which guarantees an exponential decrease for r-^oo if Im mi 
differs from zero; the function (10) then represents an approximation in accordance with the 
radiation condition at infinity. In part II of this paper (10) will appear as the first term of a 
complete expansion for the solution II Z . 

4. Position of a Turning Point for a General Profile of the Refractive Index 

The significance of a turning-point level r=r t is obvious when r t is a real quantity larger 
than the earth's radius a. We then consider the relation 

n eS (r).r- sin r (r) = W + 1)]l ' 2 ; (1 1) 

/to 

it determines a ray trajectory in a stratified troposphere with refractive index n eff) if r(r) 
represents the angle, at any level r, between the tangent to the trajectory and the radius vector 
towards the center of the earth. On the other hand, the turning-point relation m z (r z ) = can 
be written as follows in view of (9): 

„.. w . r , = M+i)r. (12) 

Therefore, r=ri characterizes the level at which the trajectory in question becomes horizontal 
( T =7r/2). In other words, the turning-point level constitutes an altitude at which this tra- 
jectory passes continuously from a rising branch into a descending branch. 

In the absence of a turning point in the physical space a<^r<^co } a zero r x of m^r) might 
occur in the interval r<^a outside the troposphere (in accordance with some extrapolation 
there of n(r)), or also at some complex value of r. The " extended W.K.B. approximation " 
(10), defined with the aid of such a nonphysical turning-point level, will still be important 
when the latter proves to be situated near the section a<>< oo of the real axis in the complex 
r-plane. 

As a matter of fact zeros r x of this type occur for the most important modes (7), since all 
these modes correspond to complex I values situated near Z defined by 

^an e M = {lo(l n +l)} l/2 . (13) 

c 

According to (12) l itself is connected with a ray trajectory tangent ing the earth's surface, r x 
then assuming the value a. Hence we expect for all eigenvalues I near l (to be discussed in 
part II) a corresponding turning point r t near a. 

The exact position of this turning point can be derived from an expansion to be obtained 
as follows. We replace (12) by the equivalent equation: 

M*(r,)=^±^=C?,say, (14) 
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in which 



M eff (r)=4-n efi (r) (15) 



denotes a modified refractive index which includes the effects of both the curvature and the 
profile of the effective refractive index n ett . We consider the following Taylor expansion of (14): 



PiaJj\ {ri a) ~ Ha* 



M? a (a)+j: ^ {ri - a y = tSg^l, (16) 



in which the dimensionless parameters 

characterize the complete tropospheric profile. In view of (14), (15), and (9) we may replace 
(16) by 

We next assume the possibility of inverting this latter relation. The inverted series can 
be represented as follows: 

r h ~a_ mf (o) A 2 4 , . /A, Aj \ „, . ( 5 A 2 A 3 5 A| A 4 \ .,. , . 

It shows how the turning point can easily be evaluated for modes with small values of mf(a) = 
^eff(^)" C?j which, however, are the only modes of practical interest. 

5. Dependence of the Extended W.K.B. Approximation on the Profile 

In view of (10) this approximation can be represented by the expression 

rf,(r)~A, { / ,(r)}1 C ^ •g'Tfo.fr)}, (19) 

which depends uniquely on the variable: 

.u«(r)=* P m l (s)ds=k a T^ {M^(s)-C\Y' 2 . (20) 

J r i J r i s 



This variable constitutes, apart from the factor i, the exponent in the asymptotic approx- 
imation of (19) for large values of |^(r)|. Therefore, u z may be interpreted as the radial con- 
tribution to the phase of the complete wave function f t (r) Pi (cos 0). 

Unfortunately, the definition (20) of u t involves the turning point r t ; hence, the expan- 
sion (18) should be substituted in order to show the complete dependence of Ui(r) on the 
profile Meff(r). Under practical circumstances, however, this dependence can also be estab- 
lished without having to resort to the turning point at all. In fact, by applying partial integra- 
tions to (20) we obtain an expansion in which the turning point disappears in all but the last 
term, such in consequence of (14). The expansion in question reads: 



u l (r)=k$L 



+iof-- Wh) ( " 1)p jrj MMs) -^' +1/2 -^Si)? j^p \ dM ^ s) 
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This expression suggests the infinite series 



222 



-W=*^Sfff-(^8) ( " 1) ' {J,ftW - 0,,rtW,, 3I32| 



d' 



{r)Y 



(21) 



dr 



In the case of convergence of this series, we can verify, differentiating term by term, the 
relation 

du t (r) k Q a 



dr 



{Ml, (r) -civ 2 ; 



this proves the correctness of (21) when taking into account the property u l {r i ) = 0. 

As an example we consider the magnetic solution in the case of the "Eckersley profile." 
The refractive index n(r) of the latter is defined by 



n 2 (r) 



\a eS \ a eB /r 2 J 



(22) 



a ett being the effective earth radius. The corresponding effective modified refractive 
index, viz 

M.„ W =^){^+(l-£)}', 

leads to the following series for (21): 

r„ m) - k ° ar ^ i t 2 2 _2_fa-a ett V+ 1 f « 2 (a) , n 2 (a)(l-ft/g e „)-61 \Mi. 

The convergence condition here amounts to: 

r 2 I flfoH n 2 (a) J I ^ 

in view of the smallness in practice of the quantities (;r—a)/a=h/a and 8Ci=Ci—n(a) this 
reduces to the following approximative condition 



2h 2a et t &C l 



a a mjS) 



<i. 



which is satisfied for all relevant values of the height h and the parameter Ci. The actual smooth 
profiles, excluding ducts and so on, deviate only slightly from the Eckersley profile; this sug- 
gests the validity of the corresponding expansion (21) in all practical cases. 




Actual profile 
Eckersley profile 



— , r = a + h 

A0 = 0'(r') -0(r") 
Figure 1. Geometry of radio ray bending. 
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6. Formulation of the Eigenvalue Problem 

It is now well established [2] that our propagation problem concerning the two media of 
the spherical earth and its surrounding stratified atmosphere can be reduced approximatively 
to a one-medium problem by introducing a proper boundary condition at the earth's surface. 
This condition is arrived at by considering the surface impedance 

y 4tt E\\ 

c H\\ 

which depends on the ratio of the horizontal components E\\ and H\\ of the electric and the 
magnetic field at the earth's surface. The factor 4x/c is such that Z represents the ratio of the 
tangential electric field at the earth's surface, and of the density of a properly chosen two- 
dimensional current distribution on this surface ; the effect of this fictitious current distribution 
is identical with that of the actual field inside the earth. 

The value of Z depends on the field solution under consideration. In the case of a flat 
boundary, and a plane incident wave arriving from the atmosphere, e.g., Z depends on the 
direction of arrival r of the incident wave. The assumption which enables us to leave out of 
consideration the medium inside the earth altogether is that Z may be replaced by its special 
value for a plane wave at grazing incidence (t = t/2). The corresponding boundary condition, 
viz 

E\\=— Z T=T/2 H\\ at r=a, 

can be worked out with the aid of (3) and (4) fdr the electric and magnetic solution respec- 
tively. It results in a relation of the following form for the scalar II describing these solutions: 

^- (rn) =- rU at r=a, (23) 

The parameter T is then given, respectively, by 

in these expressions n{a) denotes the atmospheric refractive index at the earth's surface, and 
n e the refractive index associated with a refraction from the atmosphere towards the earth. 

The boundary condition (23) can only be fulfilled for a discrete set of modes (7). Hence- 
forth the parameter I will refer to a special complex eigenvalue fixing the order of the Legendre 
function P z (cos#) that constitutes a factor of the field of such a mode. The amplitudes of 
the modes depend on the source of the field. In the case of a vertical electric or magnetic- 
dipole point source at r=b, & = the expansion in terms of these modes becomes as follows [3]: 

n= 7r^_ (2Z+1) /i(fe)/,(r)P t {cos(ir-f»)} 

C(l2b > 5 f-(rf x )1 /f(a) sin(?rZ) (24) 

The constant B represents the moment of the dipole. Each function j t (r) constitutes a 
solution of the height-gain differential equation, that has to satisfy the radiation condition 
at infinity. 
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7. Logarithmic Derivative of the Height-Gain Function at the Earth's Surface 
The boundary condition (23) can be put in the following form for the individual modes: 

r ^ (rA) l _r (25) 



t 



Hence the quantity 

r a 



I r/» J,. 



proves to be important in connection with the evaluation of both the eigenvalues /, and of the 
denominator of the amplitude factor in (24). Its approximation for small values of h/a and of 
SCi=Ci—n(a) (see the end of section 5) will be discussed now. 

We start from the extended W.K.B. approximation (19) for rf\. Its logarithmic deriva- 
tive becomes: 

^ (r/x) i k_\ <_,_„„ M i mi{u,) (20) 

Another logarithmic derivation of the expression 

<=&{ AfA(r)-0!} 1 ' I -*o{AfA(r)-0!}>« (27) 

results in: 

< M eff (r)-M:„(r) 
«i~ 2lfA(r)-0| ' 

We next introduce a definition for the effective-earth radius for any profile, viz: 

ro eff (r) 



(28) 



J r {mAr)} 



M t „(a) . . 

= M e ' ff (a)' ( - J) 



This definition involves the following value of (28) at r = a: 

(<\ „ V*M (30) 

On the other hand, the expansion (21) may be approximated for small h/a and dCi by its first 
term. Its value for r=a reduces to: 



ha efi {nl s (a)-C V* 
3 n 2 eS (a) 






A comparison of (30), (31), and of the special value of (27) at r = a then proves the vanishing 
(for the approximations under consideration) of the quantity 

£%_1^X 

6% 2< 

at r=a. Further, when evaluating the remaining term in the right-hand side of (6) for r=a, 
we may substitute, applying (27) and (31) 



«;(o)~*„(»:„(«)-ci)""~fc-f25kW l (32 ) 



•\^S7^ W '•»"•• 
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The following final approximation thus results from (26) for r—a\ 

6 \ a eB J U " W Hi)\{u x (a)} (33 > 






tfx 



8. Equation for the Eigenvalues 

According to (5) and (33) the eigenvalues are to be determined from the equation: 



"{hbC®}" (34) 



m%{u,(a)} T 
ff${«i(<*)] 



The resulting roots ifcj(a) fix in succession the corresponding values of C t and / (see (31) and 
(14)). 

The only profile parameters entering in the equations (34), (31), and (14) for the eigenvalues 
are n(a), and the effective earth's radius a eit denned by (29). Hence the conventional diffrac- 
tion theory for a homogeneous atmosphere with refractive index n(a) can be applied at once, 
provided that the actual earth's radius a is replaced by a ett . This simple result depends on 
the smallness of 8Ci=Ci—n(a) which will be proved later on in the rigorous theory of part 
II. It then appears that the eigenvalues do only depend, in the approximations under con- 
sideration, on profile properties of the refractive index near the earth's surface; these properties 
concern n(a) and the derivative n f (a), both of which determine the parameter a ett . However, 
we emphasize the assumption of a smooth profile n(r) without any discontinuity of a derivative of 
any order; such a discontinuity would involve higher-order terms in the expansion (21) for 
Ui(r) which are not negligible with respect to the first term. 

9. Approximate Expansion in Terms of Modes 

In order to obtain a final approximation replacing (24) we have to evaluate, among others, 
the derivative of (33) with respect to X. This can be facilitated with the aid of two relations, 
to be derived as follows. First we obtain from a differentiation of (31): 



du\(a) —3 
, -/■>*/ 

dX 2 



-w^^ 



The last factor can further be reduced, with the aid of the approximation Ci~l/k a of (14). 
We find: 

d(C?)_2X+l 2{X(X+1)} 1 / 2 2<7x 



dX k\a 2 Tela 2 Jc a 

On the other hand, we shall make use of the identity 



(36) 



The application of (35), (36) and (37) to the determination of the X derivative of (33) results 
in an expression which can further be simplified for X equaling an eigenvalue L In fact, we 
then have C\-=Ci^n eS (a); this latter zero-order approximation [which is equivalent to 
l~k Q a, n ett (a)] is well known from the diffraction theory for a homogeneous atmosphere. More- 
over, we may then also eliminate the Hankel functions with the aid of (34). The procedure 
outlined here results in: 






3 2/i H l3 a l J 3 n eB (a) 



b m ^+Hmm} m } (38) 



r=a;\=l 
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We further substitute in (24): 



2l+1^2l<^2k an eS (a) 



(39) 



as well as the well-known approximation for Legendre functions with large complex order 
having a positive imaginary part. The latter reads: 



P z {cOs(7T-fl)j 

sin (irl) 



<-{2i) 1 ' 



Aid 



(jl sin tf) 1 ' 2 

In view of the approximate value l^k an eU (a) we introduce the representation 

l=k an eS (a)+5l, 



(40) 



(41) 



whereas the correction 81 may be omitted in the factor Z 1/2 . The substitution of (38), (39), 
(40), and (41) into (24) leads to: 



u^- 



' cba^{an eS (a) sin tf} 1 ' 2 ^ 11(a) 



p i bid 



l><<«»"+5{j!ife}'"] 



(42) 



We arc interested first of all in the field strength attenuation that is due to the atmospheric 
refraction and the diffraction by the earth. Therefore, we next pass to the ratio of |11| and 
the modulus [n^l of the scalar determining the primary unattenuated field. The latter modulus 
reads [4]: 

in,,!-^- 1 



cb TP cba$ 

TP being the distance from the transmitter to the receiver. 
We thus obtain from (42) and (43): 



(43) 



I — 

n„ 



2 3 ' 2 k, 



* /2 h m r wad Y' 2 



s 



M)Hr) 



m«) 



[^^<AwmYl 



(44) 



The effects of the elevations b— a of the transmitter, and r— a of the receiver, are con- 
tained, for each individual mode, in the height-gain factors 

■fi(a) 1i{a) 

These factors depend, according to (19) and (20), on the complete profile between the earth's 
surface and the levels of the transmitter and receiver. On the other hand, these factors can 
be left out of consideration when considering propagation along the ground. In this case the 
field is completely determined by n ett (a), a ett and T; the dependence on the profile of the re- 
fractive index (if assumed as smooth), then only concerns its properties at the vicinity of the 
earth's surface. 



10. Connection of the Height-Gain Factors With Horizon Distances 

For an investigation of the height-gain effects we can restrict ourselves to an elevated 
receiver, the dependence on the height being completely similar for an elevated transmitter. 
Substitution of the extended W.K.B. approximation (19) in both/ z (r) and jfj(a) yields: 



/»(r) 
Ma)' 



af u^-uija) V /2 Hj)l{u t (r)} 
'r\«,(o).«{(r)/ Hi}l{u t (a)}' 



(45) 
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The effect of the height h=r—a enters implicitly through the variable Ui(r), and also in 

the factor 

1 1 

r{u , l (r)}^~(koar)^{Ml s (r)-C 2 l } 1/ ' 

The dependence on the height can therefore be expressed by: 

Mr) „ FM)} { . 

f l {a) r^{Ml s (r)-CW* { ^ } 

The variable u t (r) can be connected with a geometric property of the profile, viz the 
angular distance #(V) from the elevated receiver to its horizon point. This distance can be 
evaluated with the aid of SnelFs law for the curved ray trajectory that passes through the 
receiver, and meets the earth tangentially. For our spherically symmetric medium this law 
amounts to the relation (11), if its refractive index is given by n eil (r) instead of n(r) [thus 
differing slightly from the actual refractive index in the case of the electric solution, see (6)]. 
We have: 

r'n^ (r') sin r{r f )=an en (a), (47) 

if r(r') represents the angle, at an arbitrary point (?',#') of the trajectory, between the tangent 
of the latter and the vertical through this point. Hence tan T = r f d$\dr t '. 
The relation (47), or the equivalent one 

M eB (r') sin r(r')=M eff (a), 

can also be put in the form: 

rM / _tanr(r / )_ M eff (a) 



dr' r' r'{Ml,tr')-MUa)Y 12 



An integration along the trajectory, from the horizon point at r f = a up to the receiver at 
r'=r, yields the following formula for the angular horizon distance #(r) of the receiver: 

C r dr f 

tf(r)=M ea (a)l r>{Ml^)-M\{a)Y" <«> 

We also need the derivative of this function, viz 

d* M eU (a) 

dr r{Ml s (r)-Ml s (a)} 1/2 ' ^ } 

The connection between r and #(r) is unique for the smooth troposphere profile (excluding, 
e.g., "ducts") under consideration. Therefore, all functions of r may also be considered as 
functions of #(r), in particular the variable u t (r). In order to determine the explicit dependence 
of u t (r) on &(r), we first derive, with the aid of (49), the relation: 

d(M 2 ) 
d{M 2 eS (r) -MIM } 1/2 _ r dr 



d& ~ 2n(a) 



We can deduce similar relations for the higher-order derivatives of [M 2 eS (r)—Ml s (a)} 1/2 with 
respect to #, again using (49); for instance we find: 

| 2 {MU.)-M» B ( a )}-= 2 ^|{/if)}.(My r )-MM«)}'«. 

The special values of these derivatives at r=a can be expressed in terms of the coefficients 
Aj defined by (17). These latter derivatives determine the coefficients of the Taylor expansion 
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of [M 2 eS (r) — M* ff (a)} 1/2 , with respect to #, taken at the point r=a, that is tf(r) = 0. The expan- 
sion in question starts as follows: 

Myr)=Mya) +i J|_ ^)+f^ *<r)+ . . . ; (50) 

moreover, the definition (29) implies the relation 

A 1= 2r4(«) ?-. (51) 

We next consider the approximation of u z (r) by the first term of (21), viz 



Ui(r)^-ha 



2, {Ml e (r)-C^ 2 2 t a {MM/O-tffF'WM'O-MMa)} 1 



3 "^ r d(Ml«) 3 rc ef f (a) d(Ml s )/d# (52) 

Eacli factor can be expanded with the aid of (50), but we shall only retain terms up to the 
first one depending on &. This requires the validity of the relation 

(A 1 +A a ).f»*(r)<n a (a). (53) 

Moreover we substitute [compare (14) and (41)]: 

C?={n en (a) + 5C\}^7iUa) + 2n en (a)-8C l ~nlM+2n en (a)^ a - (54) 

The corresponding evaluation of (52) results in: 



u>i(r)' 



3{n eff (a)} 1/2 



{2 />/ a 2 "\ 3/2 

-^+n e „(a) £- ■ d*(r) V ■ (55) 

ft c,a a eS j 



Apart from the inequality (53) this approximation depends on the neglect of the higher- 
order terms of (21). This neglect proves to be justified in view of the smallness of 5l/k a and 
# 2 compared to unity, remembering also the order of magnitude of unity for n ett (a) and a/a en . 
Moreover, (53) is even satisfied automatically for all t?(r)<«l provided that we have |A 2 K<Ai- 
Under such circumstances u t (r) only depends on the eigenvalue 81, the horizon distance #(r), and 
the profile properties n cft (a) and a e ffj on the other hand, the eigenvalues 5/ merely depend on 
n ett (a) and a ett (see section 8). As a consequence n ett (a), a eit and #(r) may be considered as 
the only quantities determining u t (r). The same then holds for F{ui(r)} in (46), as well as 
for the quantity [compare (50) and (54)], 

{£/ A 2 ^ !/ 4 

The still remaining factor r~ 1/2 of (46) can be approximated by a~ 1/2 . We thus finally conclude 
that the height-gain factor does depend, in the approximations under consideration [the most 
restrictive one being (53)], on no other quantities than n efi (a), a eU and &(r). 

1 1 . Asymptotic Approximation of the Height-Gain Factors 

We might substitute the approximations of the preceding section into the height-gain 
factors (45) for both the transmit ter and the receiver, and next the values of the latter into the 
mode expansion (44). This would show the explicit dependence of the attenuation factor 
In/II^I on the four quantities n ett (a), a ett , &(r{) and #(r 2 ); #( r i) and &(r 2 ) here represent the 
angular distances from the transmitter and receiver to their horizon point. The complicated 
final expression depends, among other things, on Hankel functions of order 1/3. The latter 
can be simplified considerably by replacing them by their asymptotic approximation. This 
can be justified in the case of an argument with modulus well above unity. Unfortunately, 
the argument Ui(a) of the Hankel functions in the denominator of (45) proves to be of the 
order of unity, but the other argument u t (r) of the Hankel functions in the numerator of (45) 
may be much larger for elevated transmitters and receivers. 
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Let us consider, once again, the height-gain function for the receiver, and estimate the order 
of magnitude of \ui(r)\. When leaving out of consideration the term with # 2 in (55) we find 
an order of unity for \ui(r)\, remembering (from the case of a homogeneous atmosphere) the 
order of (k a) m for 81. Therefore, |i^(r)|>>l will certainly hold provided we have 

Hence, a sufficient condition for applying the as3miptotic expression for H[%{ui(r)} reads 
very roughly [n ett (a) and a/a etI being of the order of unity]: 

This inequality holds for altitudes h=r—a above some critical value, h cr say. 

Let us now assume K^>h cr for the receiver in question. We may then also replace the 
expression between braces in (55) by the first two terms of its binomial expansion with respect 
to U, so as to obtain: 

«,(r) ~ hn ^f a " #*(r)-5l-M& (58) 



1/2 e*{»'( r >-fH 
Hi)l{u,(a)} 



The substitution of the asymptotic expression for H[%{ui(r)} into (44) first results in: 

ji 0) ~ r \ t u t (a) -Ui (r) / 

In the exponent we must take into account both terms of (58), but the dominating first 
term suffices in the nonexponential factors. We shall substitute, moreover, the relations 
(31) and (32) for evaluating Ui(a) and u/(a) [applying also (54)], the approximation a/r^l, 
and finally the following relation based on (20), (50), (54), (51) and (56): 

^ KtfL CLeff J a ett 



ffcnW off (a)a 3 5 

i\ °f -tf 3 (r)-3/.g(r)— - 



We then find: 

jjM /ey 2 «j 

Ma) W r 2 ^ M ;, 2fJ aT fl«(-2*Z)»/' 1 ^ 

{ 2«tf(r)} ' H v*l S {k^(a)y>W*j 

12. Mode Expansion for Transmitter and (or) Receiver Above the Critical 

Altitude 

We first assume a receiver at an elevation h 2 =r 2 —a above the critical height h cr given by 
(57), and a transmitter on the ground. When evaluating the mode expansion (44) we can 
then substitute (59) for the height-gain factor of the receiver, that of the transmitter being 
unity. The term proportional to # 3 in the exponent of (59) drops out, it being real and 
independent of I; it therefore has no effect on the modulus |n/n pr |. 
The use of this substitution results in: 



II 



n pr 



2" /2 kl /e f 3afl 1 1/2 

' aU? \n e Aa)Hr 2 )J 






^( 5 0-[{3 M Ka)}-+|i{^ y } 2/3 )ff 1 'H{^(a)} 

(60) 
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in which 



«i(a)= ; 



a eff (-26l) 3 ' 2 



We next consider a transmitter and receiver at elevations hi=r\ — a and h 2 = r 2 —a which 
are both above the critical height. We may then apply (59) to each of the two height-gain 
factors fi(ri)/fi (a) said f i(r 2 ) /f i(a) occurring in (44). It leads to the expression: 



n 



ii 



gjJ.l/6 



(2a#) !/ 



1/2 X \2-J 



J8l\ z '>-»(r l )-d(r 2 )} 



^[(3^)}-+5{^} 2/3 ].[H| 



)\Mo)] 



(61) 



13. Final Expressions for the Vertical Field Component 

So far we have derived expressions for the scalar quantity II. The transition to the field 
components is conditioned by (3) and (4). In the case of the electric solution this transition 
is effected for the vertical component E T of the electric field by multiplying the coefficient of 
each individual mode by [5] il(l-\-l)/(kdr). In view of our approximations this quantity can 

be simplified to 

il 2 

By combining this quantity with the expression (43), that is 

In I B B 

1 pA eba§ ca 2 § 

we deduce the following converting factor for passing from the attenuation factor |n/n„ 7 | to \E T \: 



k a.nl n (a)-3-= WM 



C(l 2 # 



ca§ 



A similar factor proportional to # -1 applies to the reduction of |n/n pr | to the vertical magnetic 
field component \H T \ in the case of the magnetic solution. 

With the aid of this factor we find the following representation replacing (44): 



\E,\ 



2 3 < 2 kl> e >Bnl< t 2 (a)/ w 



call, 3 



w) ) 



x 



^ /,(6)/,(r) 



[i^< /3 +5{^)} 2/: 



(62) 



We also mention the field representations corresponding to (60) and (61) for a transmitter 
and (or) receiver well above the critical height: 



,„. WBiMja) f 3 \ 1/2 V 
1 r| ~ calf \a#-Hr 2 )f X 



e i 5l{#-d(r 2 )} 



(hi=0]h2^>h cr ) 
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(63) 



and 



\Er\ 



6Z#W eff % 



ca 



1/3 
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ira#-#(ri)-#(r 2 ) 
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(A 5Z { ^ — ^ (r x ) — ^ (r 2 ) } 



(64) 
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All these expressions contain the quantity H(i%{ui(a)} either explicitly, or implicitly in the 
height-gain factors. This quantity can not be reduced any further for the most general soil 
conditions of the earth. However, a convenient expression may be obtained in the two limit- 
ing cases of very long and very short waves. The latter case corresponds to an infinite value 
of the right-hand side of (34). Let us denote this right-hand side by —e iir/d N. The solution 
of (34) can then be expanded with respect to iV" 1 , a method well known in the diffraction 
theory for a homogeneous atmosphere [6]. A corresponding expansion then exists for 
H{)l{ui(a)}. Its leading term for large \N\ reads: 

tti/si u iW J~ j^ > 

u hm here marks a zero of H{)l(u). 

14. Comparison Between the Fields Corresponding to Two Different Profiles 

As stated above (end of section 9), the influence of the refractive-index profile is restricted 
to that on the parameters n eS (a) and a eff for propagation over the ground. Moreover, the 
height-gain effects reduce to an additional dependence on the angular horizon distances #( r i) 
and #(r 2 ) of the elevated transmitter and receiver (see section 10). This involves an extremely 
simple relation between the fields corresponding to two profiles having identical values of 
n eH (a) and a eff . Such two profiles show the same surface values of the refractive index and 
its gradient, but they may differ noticeably at high altitudes. According to section 8 the 
eigenvalues l ) too, prove to be identical for these profiles; therefore, the same then holds for 
all parameters occurring in expansions such as (62), (63) and (64). 

We shall compare the actual profile with some reference profile, the quantities of which 
will be marked by a dash (see fig. 1 in which the Eckersley profile (22) is assumed to be the 
reference profile); both profiles are assumed to have coincident values of a and a eff . In view 
of the above remarks, raising of an antenna (transmitter or receiver) to an elevation h=r—a 
will produce the same field as an elevated antenna at a height h'=r'—a in the reference profile, 
provided that we have &(r)=$'(r r ). 

This may be applied to atlases of propagation curves, such as edited by the C.C.I.K. [7] 
and by the Japanese Ministry of Postal Services [8]. The data of these atlases are essentially 
based on the Eckersley profile (22) with a eff =(4/3)a. The evaluation of (48) for the angular 
horizon distance results for this profile in the expression: 

1/2 



*Ec k (r)=(^) 1/2 arccos0) 



Therefore, the curves of the mentioned atlases can also be used for other profiles with the same 
value of a eS [the value of n eS (a) is irrelevant, since n eS (a)^l]. The antenna heights h=r—a, 
corresponding to a horizon distance $(r), are then to be replaced by an effective height h' —r f — a 
for which 






cos 



Examples of the use of such corrections are given in two recent publications by Norton [9, 10]. 

15. Concept of the Effective or n Angular" Distance 

The most striking feature of the expressions (63) and (64) concerns the presence of the 
angular distances #— #(r 2 ) and &— &(ri) — #(r 2 ) in the exponent. The effect of the antenna 
heights is mainly contained in the corresponding exponential, the influence of the nonexpo- 
nential factors {#( r i)} 1/2 and {#(r 2 )} 1/2 being much smaller. The importance of the distances in 
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question has been emphasized in particular by Norton, Rice, and Vogler [11]. These authors 
have introduced the term " angular distance" with an extension to the cases of propagation 
over irregular terrain. For propagation over a smooth spherical earth beyond the line of 
sight the "angular distance" represents simply the actual distance between the radio horizons 
of the transmitter and the receiver divided by the radius of the sphere. 

For a smooth earth the expansions (63) and (64) justify the reduction of t lie height depend- 
ence to that on the "angular" distance, under the following conditions: (a) the elevated 
antennas are to be situated above the critical height h cr fixed by (57); and (b) the eflfecl of l lie 
square roots {#(ri)} 1/2 and {#(r 2 )} 1/2 should be negligible. 

Let us abstract from the factor # 1/2 in (63) and (64) which is not ver\ critical. The 
dependence on the "angular distance" then combines both the influence of the antenna eleva- 
tions and of the horizontal distance. The influence of the former simply amounts to a reduction 
of the actual distance d=a$ to that of the radio horizons, that is to a{& — #(/'i)} or a{&— #( r i)"~ 
#(r 2 )}. The comparison with a reference profile (see the preceding section) can be expressed 
here in terms of a shift M of the "angular" distance. In fact, lei an antenna height h=r—a 
in the actual profile produce the same field, at a distance d=a&, as the reference profile at a 
distance d-\- A =a(# + A#), the antenna elevation being kept constant. The unique dependence 
on the "angular" distance, assumed here for both profiles, then involves tin 4 relation: 

#— #(r)=#+A#— #'(*•). 
Hence 

Atf=tf'(r)— 0(r). 

The corresponding linear distance shift is given, according to (48), by 

r r dr' r 1 1 "I 

A=aAtf =aM eff (a)J a y- [ { M: ff »(rO-Mi(a)} 1 / a "{M! ir (r')-Mi r (a)}^J 

The derivative with respect to the antenna elevation h=r—a, reading 



dA 
dh 



-MJfi) [ { M3 (r) _ ML {a)}m {MIb {r) _ M ^ (o) ]m j 



can be reduced for small altitudes with the aid of the expansions: 

Ml R (r)-M! n (a)==^h+^+ . . ., 

M:£{r)-MlM=^h+!^K+ . . . 

This leads to the following approximation for A itself for small heights: 

n eff (q)(A 2 -A;) /2 (A 2 — Ag)q^ 2 , 3/2 
6a 1/2 A! 3/2 \2.2 ll2 nlM^ 

Equivalent expressions have been derived by Millington [12] and Wait [13] in a simplified 
derivation of the effective-distance concept. 
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